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Abstract 

We study statistical properties of zeros of random polynomials and random analytic functions 
associated with the pseudoeuclidean group of symmetries 517(1, 1), by utilizing both analytical 
and numerical techniques. We first show that zeros of the 5?7(1, 1) random polynomial of degree 
A'' are concentrated in a narrow annulus of the order of A'^"^ around the unit circle on the com- 
plex plane, and we find an explicit formula for the scaled density of the zeros distribution along 
the radius in the limit N — > oo. Our res ults a re supported through various n umeric al simula- 



tions. We then extend results of Hannay [Han] and Bleher, Shiffman, Zelditch BSZ2] to derive 
different formulae for correlations between zeros of the SU{1, 1) random analytic functions, by 
applying the generalized Kac-Rice formula. We express the correlation functions in terms of 
some Gaussian integrals, which can be evaluated combinatorially as a finite sum over Feyn- 
man diagrams or as a supersymmetric integral. Due to the 517(1, 1) symmetry, the correlation 
functions depend only on the hyperbolic distances between the points on the unit disk, and we 
obtain an explicit formula for the two point correlation function. It displays quadratic repulsion 
at small distances and fast decay of correlations at infinity. In an appendix to the paper we 
evaluate correlations between the outer zeros \zj\ > 1 of the 517(1, 1) random polynomial, and 
we prove that the inner and outer zeros are independent in the limit when the degree of the 
polynomial goes to infinity. 



1 Introduction 



In this paper we are interested in statistical properties of zeros of random polynomials and random 
analytic functions associated with the pseudoeuclidean group of symmetries 5L/(1,1). The moti- 
vation for the study of zeros of random polynomials and random analytic functions comes from 
different applications, most importantly from the theory of quantum chaos (see papers of Bogo- 
molny, Bohigas, and Leboeuf |BBL|, Leboeuf and Shukla |LS], Hannay |HaE], Korsch, Miller, and 
Wieschcr [KMW|, Nonnenmacher and Voros |NV|, Forrester and Honner [FH|, Leboeuf Q], Shiffman 
and Zelditch |SZ1|, and others). There are different ensembles of random polynomials associated 
with different groups of symmetries. In this respect we can make a parallel to the theory of random 
matrices. In the latter, the basic ensembles include orthogonal, unitary, symplectic, circular, and 
some others (see |m|). Making the parallel, we may think of the 0(n -I- 1) ensemble of real poly- 
nomials as of an analogue of the orthogonal ensemble of random matrices. The 0{n + 1) ensemble 
of random polynomials consists of multivariate homogeneous real polynomials of (n + I) variable 
z = (zo, zi, . . . , z„) of the form 



^(z) 



(1) 



\m\=N 



where m = (toq, it^i, ■ ■ • j "m-n) is a multiindex, \m\ = niQ + mi + • • • + to„, z™ = z™'^ . . . z^", Cjy 
is the muhinomial coefficient, 



C 



^ TOo!mi! . . . m„! 



(2) 



and ttm are independent standard, N(Q, 1), real Gaussian random variables. Consider a set Z of 
joint real zeros oik, k < n, independent copies of the 0{n + 1) random polynomial in the projective 
space RP". This is a random real algebraic variety of dimension n — k. It is nondegenerate almost 
surely and it possesses a natural volume element induced by the standard metric in MP". The joint 
distribution functions of the zeros are invariant with respect to the action of the group 0{n + 1), 
and, in particular, the density function of the zeros is constant (cf. phSm |, [BSZ2], [BSZ3|). 

Similarly, the SU{n + 1) ensemble of random polynomials consists of multivariate homogeneous 
complex polynomials of (n + 1) complex variable z = {zq, zi, . . . , Zn) of form (p]) w here a.^ are 
independent standard complex Gaussian random variables (cf. [BBL|, |Han|, [BSZ1]-|BSZ4|). It 
corresponds to the unitary ensemble of random matrices. The distribution of joint zeros of k in- 
dependent copies of the SU{n +1) random polynomial is invariant with respect to the action of 
the group SU{n + 1) on the complex projective space CP". As the degree N of the SU{n + 1) 
random polynomial goes to infinity, the scaled correlation functions of zeros approach a limit, which 
i s rep r esented by the correlation functions of the Wn ensemble of random analytic functions (see 
[ Han| , BSZ1 |-| BSZ4 |, The latter consists of random functions of the form 




(3) 



where the sum runs over multiindices m = (mi,...,m„) with rrij > 0, ml = mi!...,TO„!, and 
are independent standard complex Gaussian random variables. The existence of the limit for 
the scaled correlation functions of zeros is valid in a very general framework of random sections 
of powers of line bundles over compact manifolds and the limit is universal (see [BSZl|-[ i3SZ4[ ). 
The explicit combinatorial formulae for the limit of two-point correlation functions are obtained in 
[ BSZ2| -|BSZ4|. They demonstrate a quadratic repulsion if fc = n = 1, neutrality if fc = n = 2, and 
attraction if fc = ri > 2. 

In this paper we are interested in the pseudoeuclidean, SU (1,1) ensemble of random analytic 
functions. The general SU(n, 1) ensemble consists of multivariate analytic functions of the form 



i^iz) 



E 



^[fa\+L-l "i 



(4) 



where L > 1 is a fixed integer, a parameter of the ensemble, the sum runs over multiindices m 
(toi, . . . , m„) with mj > 0, = zf^ . . . z™-, 



C" 



i\ + L-l 



{\m\ + L-l)\ 
[L — 1)!toi! . . . m 



(5) 



and a„i are independent standard complex Gaussian random variables. We will restrict our study to 
the case n — 1. The basic calculations are extended to the case of any n, and we are going to return 
to this extension in subsequent publications. Another possible extension concerns the universality 
result like in |BSZ2| , [|z|, and jBSZaj . 

The plan of the paper is the following. In Section || we consider the ensemble of SU{1, 1) random 
polynomials, which is obtained by restricting m in (^ to be bounded by N (with n = 1). We calculate 
the density function for the distribution of zeros of the SU{1, 1) random polynomial and we show 
that the zeros are concentrated in a narrow annulus around the unit circle, of the width of the order 
of 1/N . We find the scaled profile of the density function along the radius in the limit N oo. In 
Subsection 2.2 we rederive the result of Leboeuf Q for the density of the SU{\, 1) random analytic 
function. In Section |^ we derive different formulae for the correlation functions between zeros of the 
5?7(1, 1) random analytic function. First, we apply the general result of [ BSZ2| to get the Kac-Rice 
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type expression for the n-point correlation function. We then specify it for the case of a Gaussian 
random field (cf. |Han|, pD( |, |BSZ2|). This expresses the n-point correlation function in terms 
of some Gaussian averages, which is further represented as a sum over Feynman diagrams, or in a 
different approach developed in |BSZ4|, as a supersymmetric integral. In Subsection 3.3 we carry 



out concrete calculations for the two point correlation function in the SU{1,1) ensemble. They 
demonstrate a quadratic repulsion at small distances and fast decay at infinity. In the limit when 
the main parameter L of the SU (1,1) polynomial goes to infinity we recover the Wi correlation 
function of Hannay. There is an appendix at the end of the paper, in which we consider correlations 
between outer zeros Zj, \zj\ > 1, of the SU{1,1) random polynomial. We show that in the limit 
— > oo, the outer zeros are independent of the inner zeros \zj\ < 1, and after changing the variable 
z to 1/z, the correlations between outer zeros coincide with the correlations between zeros in the 
SU{1, 1) ensemble with L = 1. 



2 Basic Statistical Properties 
2.1 Scaled Density Function 

Consider the following random polynomial, associated with the pseudoeuclidean group SU{1, 1): 

where L is a fixed positive integer, a parameter of the problem, C™_|_^_^ are Newton's binomial 
coefficients, and am are independent standard complex Gaussian random variables, so that 

Earn = 0, EamOn = <^mn, EamO-n — 0. (7) 

When L = 1, (^) reduces to the classical form 

N 



m=0 



We wish to investigate the density of zeros pn[z) for the SU{1, 1) polynomial (||) as iV ^ oo. The 
density pn{z) is determined by the condition that for any test function ^p^z), which is infinitely 
differentiable and compactly supported, 



E 2.f{zj) = / PN{z)f{z) dz, dz = dxdy, 



(9) 



where Zj are zeros of the random polynomial ^{z). Since the total number of zeros is equal to N, 
we have that 

/ PN{z)dz = N. (10) 
Jc 

Observe that the polynomial 

N 



m=0 



i;{e''z) = J2 V^™+L-i ame"^'^"" (H) 



has the same probability distribution as ip{z), hence 

PN{r^')^PN{r). (12) 
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The general formula for pn{z) is given by the Poincare-Lelong type expression 

1 d'^ 



Pn{z) 



TT dzdz 



[In E{^{z)ij{z))] 



(see PSZII , U). From (g), 



N N 
m— 7n— 

In particular, 

AT 



where here and in what follows we use the notation 



1 ^ 

— y E (™ + 1) • ■ • (™ + ^ - ^W' 



m=0 



X = zz = \z\ 



It is obvious that 



^n,l{x) > 0, Va; > 



Substituting ( ]l5| ) into ( p^ ) gives that 



We further define 



PN^z) 



N 



Then from (|T|), 



Qn,l{^) = E 



r;^+l - 1 
X - 1 



1 c;^~'gjv,L(x) 

{L-l)\ dx-^-i 



(13) 
(14) 

(15) 

(16) 
(17) 

(18) 
(19) 
(20) 



Numerical simulations show that most of the roots of SU{\, 1) polynomials (g) are concentrated in 
a small annulus (of the width of the order of 1/A^) near the unit circle. This property is illustrated 
in Fig. 1, which contains a point-plot of the zeros of 150 SU{\, 1) polynomials of degree N = 200, 
with L = 30. 



To get the asymptotics of the density pn{z) near the unit circle, we introduce a scaling of the 
variable x in the form 



Then, 



x = l + ^. (21) 



^=N^. (22) 
ax as 



In subsequent calculations, we will assume that 

-A<s<A, (23) 

for some arbitrary fixed A > 0. The notation R{s) = 0{N^^) used below means that there exists 
some C{A) > so that \R{s)\ < C{A)N^'^ for all s e [-A, A]. The main result of this section is the 
following theorem. 
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L = 30 
N=200 



|z| = 1 





S(z) 



Figure 1: Zeros of SU{1, 1) Random Polynomials 



Theorem 2.1. As N ^ oo, 



+ 0(iV-i) . 



where 



dsi 



(24) 



(25) 



Remark. Observe that one N in the normalization of pn{z) on the left is due to the rescaling 
( pT| ) and one is due to the integral condition (|o|). 

Proof. By substituting (^ll) into (|l9|) we obtain that as — > oo, 



1 (i + ^)A^+i_i 



{1 + 0{N-')) ^Ng{s){l + OiN-')) , 



N 



where we define 



Observe that for j = 0, 1, 2, . 



p^o + ^' 



— CX) < S < CX) 



It is obvious from (M) for s > 0. For negative s, use the identity 



1 - e-" 1 



1! 



(26) 
(27) 

(28) 
(29) 



Formula ( P6D holds obviously in a small complex neighborhood of the segment —A < s < A. 
Therefore, we can differentiate it in s, so that 



^!^^^iv^+V^)(.)(i + o(iv-i)) , 



(30) 
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Thus, from (20), we obtain that 



^n,l{x) 



{L-l)\ 



and 



hence 



(i + o(iv-i)) 



(31) 



(32) 



(33) 



Let us go back to formula (|T^. In this formula, x = 1 + 0{N~^), and the first term in the brackets 
has an extra derivative so this is the leading term, and wc can neglect the second term. This gives 
that 



1 



5(i-i)(s) 



+ 0{N~^), zz = l + 



N 



which was stated. Theorem 2.1 is proved. 



Equation ( |24D can also be expressed in the following form 

1 

TT 

From 



AT- 



N-^Pn{z) = - logg(^-i)(s)J +0{N-^) 

gW(s) 



hm Ar-p^((i + _)e-) = p{s) ^ - 



2N 



5(i-i)(s) 



log5(^-i)(s) 



Consider now the distribution function of zeros. Define 

PN{x)=N-'Eif{j:\zj\^ <x}, 
which gives the expected value of the fraction of zeros in the disk of radius y/x. Then 

P^{x) ^ N-\pn{z), x = |z|2. 



P, 



N 



Hence (24) implies that 

As follows from (p9|), the limiting distribution function, 

lim Pat f 1 + 41 = P{s 



+ 0{N-^). 



has the following asymptotics at — oo: 



N 



P{s) = -- + 0({-s)'^-^e' 
s ^ 



s — oo. 



At +0O we use the formula 

gW(.) = ^ 
s 



si 



si+i 



(34) 



(35) 

(36) 

(37) 
(38) 

(39) 

(40) 
(41) 
(42) 
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which gives that 



, , 1 L - 1 
Pi-^) = — 



For L = I, ( [40D reduces to the well-known result (see |BS|) 

^ ' sie'' - 1) 

For L = I, the limiting distribution is symmetric, so that 

P{1 -s) = I -Pis) 



(43) 



(44) 



(45) 



(it is related to the symmetry of the zeros oiip{z) in (||) with respect to the transformation z I/2). 
Since 



^(^)-^+2! + 3! + 4! 



we have that 



Thus, 



j + 1 l!(j- + 2) 2l{j + 3) 3!(j+4) 



1 



(46) 



(47) 



L + 1 l!(L + 2) 2!(i + 3) 



0{N-') 



(48) 



L 1!(L + 1) 2!(L + 2) 



In particular. 



L + 1 



+ o(iv-i). 



(49) 



which means that the expected value of the fraction of zeros inside the unit disk is asymptotically 
equal to L/ (L + 1). 

In our numerical simulations we generated a large number of SU{1, 1) random polynomials of 
degree N, calculated the zeros using standard techniques, and counted the number of zeros in annuli 
of fixed width, which are concentrically spread around the origin, per one polynomial. The scaled 
density was obtained by dividing each of these numbers by the area of the corresponding annulus, as 
well as by adjusting it with respect to N (according to the previously established result). Figure 2 
shows the results of this procedure for L = 4 and N = 150, in comparison with the theoretical value 
as given by Eq. (p^). 

2.2 SU{1,1) Ensemble of Random Analytic Functions 

In this section we consider the SU{1, 1) random analytic function. 



m—O 



(50) 



which is obtained from (g) by setting N = 00. Here again C™_|_^_]^ are Newton's binomial coefficients, 
and independent standard complex Gaussian random variables. 
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0.9 0.92 0.94 0.96 0.98 1 1.02 1.04 1.06 1.08 1.1 



Figure 2: The Scaled Density Function - Theoretical Limit and Computer Simulation 



Proposition 2.2. Random series ^5(\ ) converges almost surely for all \z\ < 1. For all mu 
distinct zi, . . . , z„, in the disk |z| < 1 the covariance matrix 

An = (e U{Zp)^{Zjy))) 

\ \ / / p,p — l,...,n 

is positive definite. 

Proof. Almost sure convergence. Observe that (m + j)/m < j + 1 hence 

(m + 1) . . . (to + L - 1) 



< LI 



and 



_ (m + £-l)...(TO+l) .L-l 
'^m+L-1 - (L-iy. - ■ 

Consider the set 

An — {a — (ao, ai, . . . ) : \am\ < m Vto > N}. 
Then for all a G series ( pO| ) converges and 

lim ProbAAT = 1, 

JV— »CX3 

which implies the almost sure convergence. 

Positive definiteness. Consider the quadratic form 



An{tj) ^ ^ E (^'ip{zp)ilj{zp')^ Hp^j.p> = E 
p,p'=i 

Assume that An is not positive definite. Then, for some nonzero vector — (/ii, . . . , 

n 

^fipip{zp) = 



for almost all ■0. Hence, ( pTj ) holds for almost V' with coefficients vector a e Aat. Assume that 
N > n. There exists a polynomial 



m=0 



such that 



so that 



(58) 

(59) 
(60) 



p— 1 p— 1 

Let a° = (flo, a?, ■ • ■ , a°_i, 0, 0, . . . ). From (|5^), it is obvious that if a G Ajv then (a + a°) S Aat. 
The shift a — > a + a*^ is a measurable one-to-one transformation in Ajv , hence for almost all vectors 
a G Ajv we have that (|5^) holds along with 



But this implies that 



which is in contradiction with 
sition 2.2 is proved. 



^^ip[il){zp) + V'o(2p)] = 0. 



p=i 



(61) 



(62) 



This contradiction proves the positive definiteness of ■ Propo- 



Remark. The above proof gives also the positive definiteness of the covariance matrix An corre- 
sponding to the random polynomial (||), provided N > n. 

In this section, we will be interested in the density function p{z) of the distribution of zeros of 
(0) in the disk |z| < 1. We have from M) that 



E{^{z)i:{z')) = C"+L-i(^^')^ 



m=0 



In particular, 



E{^j{z)i,{z)) = fix) 



(1 - zz')^ 



X = zz 



[l-x)^'' 

In terms of / the density function p{z) of the zeros of ([50|) is given by: 



Substituting ( p^ ) into this formula gives that 

L 



p{z) 



l-x 



L 



l~x 



L 



7r(l-x)2 ' 



or going back to the complex variable z, we obtain the result of Leboeuf [Q, 

L 



(63) 



(64) 



(65) 



(66) 



(67) 



Numerical simulations were similar to those of the previous section. The density was simply obtained 
by dividing the number of counted zeros by the area of the corresponding annulus. Figure 3 shows 
the results of this procedure for the degree N equal to 50, 75, 100, and 150, in comparison with the 
theoretical value as given by Eq. (pT 
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A/= 50, 75, 100. 150 
(from bottom to top) 

p(r) = (L / 7r)(1 - r^) - 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

r 

Figure 3: The Unsealed Density Function - Theoretical Limit and Computer Simulation 

3 Correlations between Zeros 

3.1 The Symmetry 

We have the following important property: 

Theorem 3.1. The distribution of zeros of the random analytic function is invariant with 
respect to the action of the group SU{1, 1), 

az + b 



(68) 



Remark. This implies that all joint distribution functions of zeros of ip{z) are ^[/(l, 1) invariant. 
Proof. Consider the following homogeneous analytic function of two variables: 



^{zo,z,)=z-^J2y^'^ 



m=0 

Then 

Let us find the covariance function of ^'(zo, zi): 

-L °° 



l^il < l^ol 



= (zoZq) 



ZlZ-^ 



= [Zqz'q - Ziz[ 



(69) 

(70) 

(71) 
(72) 



The action of a matrix A G SU{1, 1) preserves the (1,1) indefinite scalar product zqz'q — ziz[. Thus, 
the covariance function of the Gaussian random analytic function "^{zq, zi) is SU{1, 1) invariant. In 
addition, 



£;*(zo, zi) = 0, E-^izo, zi)«'(4, z[) = 0. 



(73) 
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For the Gaussian random function, the first two moments determine it uniquely. This imphes that 
the distribution of zeros of '^{zq, zi) is SU{1, 1) invariant. Restricting this to zo = 1 we get that the 
distribution of zeros of ip{z) is SU{1, 1) invariant as weU. Theorem 3.1 is proved. 
Recall that a general formula for a matrix A G SU{1, 1) is the following: 



A 



a b 
b a 



\b\'-\a\' = l, 



so that it depends on three real parameters. An SU{1, l)-invariant metric is 

(I — — y^)^ ' 

and the corresponding distance r(zi, Z2) on the disk {|z| < 1} is determined by the equation 

tanh(r/2)= 

|1 - ZiZ2\ 

The corresponding SU{1, l)-invariant volume element is 

Adxdy 
(1 — — 2/2)2 



(74) 



(75) 



(76) 



(77) 



Theorem 3.1 implies the SU{1, 1) invariance for the normalized correlation functions. 



3.2 Correlation Functions - Preliminaries 



The n-point correlation function Kn{zi^ . . . , Zn) is determined by the condition that for any test 
functions (pi(z), . . . ,(^„(z), which are infinitely differentiable and compactly supported in the disk 
{|z| < 1}, such that their supports do not intersect, 



n / 00 \ r ^ 

n X! '^P(O) = / Knizi, Zn) Y[[(p{Zp) dz 

p=l \j=l J "'C" 



(78) 



where are zeros of the SU [1,1) random analytical function i!{z) (we denote them by C,j to 
distinguish them from the variables Zj). The sum over j is, in fact, finite because (pp{z) has a 
compact support. The general formula for _ft'„(zi, . . . , z„) is given by the Kac-Rice expression 



Kn{z) 



d(D.nio,C;z)Y[i^pCp), 

p=l 



Z — (Zi, . . . , Zn), ^ — (^1, . . . , ^„), 



where £'„(a;,^;z), x,^ € C", is the distribution density of the two random vectors 
X - (^(zi), . . . , V(z„)) , S = (Vy(zi), . . . , ^\z^)) . 



(79) 



(80) 



Formula ^ is derived in [|BSZ2[| in a much more general situation of sections of powers of a lin e 
bundle over a complex manifold, as a generalization of the original formula by Kac and Rice [Ka|, 
[H (see also [Ra^ and pD|). 



For the Gaussian random functions, formula ( |79| ) can be specialized as follows (see [Han], [BD] 
and ]BSZ2|): 



Kniz) 



1 



tt" det Ar, 



\p=i 



(81) 
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where 

A„ = (E^/j{zp)^{zp,)) , (82) 

\ / p.p' — l....,n 

and (•)a„ stands for averaging with respect to Gaussian complex random variables ^i, . . . ,^n, such 
that 

(^^f?^)p,p'=i,...,„ = ^"' E^p = 0, E^p^p,=0, p,p' = l,...,n, (83) 

where 

A„ = C„ - B*„A-'B,, (84) 

and 

Br,= (E^{zp)^7{^) , Cn= (e^'{zp)WW)) ■ (85) 

V / p,p' — l,...,n V / p,p' — l,...,n 

For n ~ I, formula (^l]) reduces to 



where 



^ . ^ 1 A^iz)Cniz) ~ B,,{z)B^iz) 



dAniz) ^ , , d'A.^{z) 



A„{z) = E^{z)^{z), Bn{z) = ^^, C„(z) = — (87) 

oz ozoz 

The 1-point correlation function Ki{z) is nothing else than the density function of the distribution 
of zeros and formula ( |86| ) is equivalent to the Poincare-Lelong type formula (|l3|). 
The normalized correlation function fc„(zi, . . . , z„) is defined as 

Kn{zi, . . . , Z„) 



Ki{z-{) ■ ■ ■Ki{Zn) 
It satisfies the following theorem. 

Theorem 3.2. The Junction fc„(zi, . . . , z„) is invariant with respect to the action of the group 
SU{l,l), 

k^{Az,,...,Az„) = k„{z,,...,z^) , yAeSU{l,l), (89) 

where 

az + b , ^ 

Az^^— (90) 
cz + a 

Remark. This implies that fc„(zi, . . . , z„) is a function of pairwise distances T{zp, Zq) defined by 

®- crn . . . 

Proof. By Theorem 3T, the distribution Kn{zi, . . . , Zn)dzi . . . dz„ is SU{1, 1) invariant. Also, 

the distribution Ki{zi) . . . Ki{zn)dzi . . . dz„ is SU (1,1) invariant. Hence their quotient (the Radon - 
Nikodim derivative), 

Kn{zi,...,Zn)dzi...dZn , , . . 

— — — — — — = fc„ zi, . . . ,z„ (91 

Ki[Zi) . . . Ki[Zn)dzi ...dZn 

is SU{1,1) invariant, which was stated. Theorem |3.2| is proved. 

The normalized correlation function can be expressed as a supersymmetric (Berezin) integral, 



■ • ■ ' ^"^ = d^ y det[/ + A„(z)^^] """^ ^''^ 

where 51 is n x n matrix, 

= (W^p)p,p'=i,...,„' (93) 

and the ?7p, fjp are anti-commuting (fermionic) variables, with dr] = Yip'^^p'^^p- '^^^ integral in ( |9^ ) 
is a Berezin integral, which is evaluated by simply taking the coefficient of the top degree form of 
the integrand dct[/+A (z)n] ■ ^"-"^ ^ derivation of formula (|9^) , sec | BSZ4 1 . 
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3.3 The Two-Point Correlation Function 



From (pOD, (|67|), (p2D, and (|85|), we can directly obtain the following expressions for An, Bn, and C„ 
in the case of a two-point correlation function ^2(^1, 



B2 



C2 = 



(l-ZrZ2)-^ (l-|z2n-^ y 

Lzi(l-|zi|2)-i-i Lzi(l-zizl)-^-i \ 

iZ2(l-^^2)-^-^ Lz2(l-k2n-^-V 

L(l + - L(l + - 

i(l + izr^2)(l - ^22)-^-" i(l + ik2p)(l - 1^21')-^-' 



(94) 



By Theorem B.2, ^2(^:1,22) is SU{1,1) invariant, that is 

k2{AzuAz2) = fc2(2i,z2), VA e 5C/(1,1). 



(95) 



The action of SU{1, 1) is transitive, hence we can move zi to the origin, and by rotation, we can 
then move Z2 to the positive half-axis. Therefore, we will assume that zi — and Z2 = r > 0. In 
that case, equations (U) become 



A. = 



Bo = 



C2 



From (|96|) and (84), we further obtain 



1 (l-r2)-^ ) 
\ 

L L 

L - 



(96) 



^iif^ L(l + Lr2)(l-.^)-^--LViB^ 
For the two-point case, formula ( |8l| ) reduces to the following expression: 

1 



(97) 



^^^2 (21,22) 



7r2 det yl; 



(A2,llA2,22 + A2,12A2^2l) 



(98) 



which, when combined with (p8[), (86), and (|67|), results in the following formula for the normalized 
two-point correlation function: 



^2(21,22) 



"(1 - r2)'''+' + ((l2 - 2L - 2)r4 + (4L + 4)r2 - l) (l - r^)''' 
+ ((L + Ifr^ - (4L + 2)r2 - l) (l - r^)^ + l] ^ (l - (l - r^)^)' , (99) 
where according to ([76|), 

.^tanh^^i^l^. (100) 
127 l-2izi ^ ' 
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Figure 4: The Two-Point Correlation Function 



When L = 1, formula (^9|) simplifies to 

k2izi,Z2)^r^i2-r^), L = 1. (101) 

Plots of /c2 (tanh (0) for L = 1, 5 and 50 are shown in Fig. 4. 

As it can be seen from the plot, the two-point correlation function goes to as t ^ (or, in 
other words, as r ^ 0). The limiting behavior can be obtained through a series expansion of (|99|). 
The following expression was obtained using Maple"""^ : 

' 2 L 4 L 36 L 72 L ^ V y , V 7 

which shows dominating quadratic behavior in the neighborhood of r = 0. Thus, there is a quadratic 
repulsion between zeros. 

We are also interested in the asymptotic behavior of the correlation function as L oo. For this 
purpose, we introduce the scaling 

r=-^. (103) 

Substituting this into (^9|) and taking the limit L —> oo, we obtain the following expression: 

,^ e^"' + (u'' -4u2 - l)e2«' + (u"* + 4u2 - l)e''' + 1 , , 

k2{u) = , (104) 



which can also be written as 



(sinh^ t + t'^) cosh t ~ 2t sinh t 

sinh^ t ' ^ 2 



k2{u) - ^ Irr^ , t = - , (105) 



and agrees with the result obtained by Hannay |HaE]. A plot of (105) is shown in Fig. 5. 

In fact, under the scaling ( |103| ), the SU{\, 1) random analytic function converges, as L ^ oo, to 
the Wi random analytic function (cf. Q]). Indeed, 



^ ^ m— 
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Figure 5: Asymptotics of the Two-Point Correlation Function as L ^ oo 
As L — > oo, the expression under the radical approaches 1/m!, so that ^{u/^fL) approaches 

oc i ^ 

V'o(") ^y^\ —, ^mu"'^ (107) 

m=0 

which is the Wi random analytic function. 

A Appendix. Correlations between Inner and Outer Zeros 

In this appendix, we consider the limit of correlations between zeros of the random polynomial (^) 
as iV oo. In the open disk {\z\ < 1}, the random polynomial (||) approaches the SU{1, 1) random 
analytic function (|5^), and the correlations between zeros of the random polynomial approach the 
ones of the random analytic function. However, according to formula (p9|), there is a 1/(L + 1) 
fraction of the zeros outside of the unit disk. The limiting correlations between those, outer zeros, 
and between inner and outer zeros are described by the following theorem. Let k^{zi, . . . , z„) denote 
the normalized n-point correlation function ( |88| ) corresponding to the parameter L. Indication of 
L is important for the theorem. Let furthermore /c^^(zi, . . . , z„) denote the normalized n-point 
correlation function for the zeros of polynomial (^). 

Theorem A.l. Assume that |zi|, . . . , \zm\ < 1 and \z.m+i\, ■ ■ ■ , \zn\ > 1. Then 

lim fc^jyl^lj • ■ • I ^n) = ^m(^li ■ • • I ^ni)^n-m(^m+l' ■ ■ • ' (10^) 
N—^oo 

SO that in the limit N oo, the inner and outer zeros become independent, and the limiting corre- 
lations between the outer zeros coincide, after the change of variable z ^ 1/ z in the argument, with 
the correlations between zeros of the SU{1, 1) random analytic function with the parameter L = I. 

Proof. We will first find correlations between the outer zeros and then prove the independence 
of the inner and outer zeros. 

Correlations between outer zeros. Consider the random polynomial (j^) and define another random 
polynomial, 

ifiiz) = , ^ z^^(z-i). (109) 



\I^N+L-l 
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Observe that if zj, \zj\ > 1, is a zero of 4'{z) then Zj ^ is a zero of (p{z) and \zj ^\ < 1. Consider ip{z) 
in the disk \z\ < 1. From (|6|), 

(f{z) = flAT + 



N ...{N + L-2) 
(N + l)... {N + L-1) 



(Ar + l)...(Ar + i- 1) ^ ^ ^ ' 



As A'^ — > oo, the expressions under the radical approaches 1 from below. In addition, we can replace 
aj^-m by a,„, because they are the same standard random variables. Thus, as iV ^ oo, ip{z) 
approaches the random function 



ip{z) = ao + aiz + a2Z^ + . . . , \z\ < 1, 



(111) 



which is the SU{1, 1) random analytic function with the parameter L = 1. Observe that ip{z) is a 
Gaussian random polynomial and for its correlations we have formula (|8^). From this formula, we 
obtain the convergence of correlations between zeros of ip{z) as A ^ cx) to the ones of the SU{1, 1) 
random analytic function with L = I. 

Independence. We introduce the random function ri{z) such that i]{z) — ip{z) for \z\ < 1 and 



1 



z-^V(^), 1^1 > 1. 



(112) 



Then the zeros of r][z) and ipi^) coincide. In addition, r](z) is a Gaussian random field and its 
covariance function E{r]{z)rj{z')) coincides with E{%p{z)il}{z')) when |z|, \z'\ < 1, and hence, as A^ — > 
oo, it approaches the correlation function (^) of the SU (1,1) random analytic function with the 
parameter L. Similarly, as we saw above in this appendix, when |z|, \z'\ > 1, E{7]{z)rj{z')) approaches 
the SU{1, 1) covariance function with L = 1, if we replace z, z' by z^^, z' ^, respectively. Consider 
now the correlation function E{'q{z)r]{z')) when |z| > 1 and |z'| < 1. From ( |112| ), 

N 



EivizUz')) 



.,N 



^l!+m-l{zz'Y 



(113) 



Assume that |zz'| < 1. Then we can estimate the correlation function as follows: 

N 



EMzHz')) 



< 



Ho 



L+m-1 



< 



z 



N 



(114) 



If |zz'| > 1, then we similarly obtain that 



EivizMz')) 



,N 



N 



'~^T+m~l{zz') 



m=0 



< A^-^k I 



Combining the two cases we can write 



E{r^{zMz')) < A^(max{|z|-Mz'|}) 



N 



(115) 



(116) 



This shows that the values of ri(z) inside and outside of the unit disk become independent as A^ — ^- oo. 
Hence their zeros become independent. Explicit estimates for the correlations between the inner and 
outer zeros follow from formula (^Tj) applied to 77(z). 
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